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Voltage-driven spin transfer torque in a magnetic tunnel junction comprising magnetic insulating 
electrodes is studied theoretically. In contrast with the conventional magnetic tunnel junctions 
comprising transition metal ferromagnets, the spin transfer torque presents unconventional bias 
dependencies, related to the presence of spin-dependent Fowler-Nordheim tunneling processes. In 
particular, we find that (i) the out-of-plane torque generally dominates the in-plane torque, (ii) out- 
of-plane torque and in-plane torque are symmetric and asymmetric at low bias voltage, respectively, 
and (iii) both of torques show a dramatic enhancement at large bias voltage. Materials consideration 
are discussed and we show that due to the low damping parameter expected in magnetic insulators 
a spin transfer torque can be experimentally observed in such systems. 



PACS numbers: 72.25.-b,73.43. Jn,73.40.Rw,73.43.Qt 

I. INTRODUCTION 

Ferromagnetic and ferrimagnetic insulators have 
gained increasing interest in the recent years with the 
observation of thermally-induced pure spin currents at- 
tributed to the so-called spin Seebeck effecP. These 
magnetic materials possess a very low magnetic damp- 
ing parameter (a ss6.710~ 5 for Y3Fe50i2) resulting from 
the absence of conduction electrons and constitute ideal 
tools to investigate spin waves propagation and non- 
linear properties over long distances^. Indeed, the elec- 
trical generation of spin-waves in YaFesOii^H (YIG) 
through spin transfer torque 5 -' as well as the obs erva tion 
of spin pumping^ from YIG into Au and Pt layerP^have 
proven that these complex materials may be relevant for 
spin transfer torque and spin pumping studies. Whereas 
Slonczewski proposed that magnons thermally excited in 
a magnetic insulator can exert a torque on an adjacent 
free layeiP^, Xiao et alP^ were able to interpret spin See- 
beck experiments^ in terms of spin pumping induced by 
thermal magnons. 

Besides their low magnetic damping constant, mag- 
netic insulators are also expected to make a difference 
in the race for the largest tunneling magnetoresistance 
(TMR) effect in magnetic tunnel junctions (MTJs) com- 
posed of two ferromagnets separated by an insulatoi^. 
While pioneering studies used amorphous AI2O3 (alu- 
mina) barriers and transition metal electrodes^, an out- 
standing breakthrough was achieved ten years ago with 
the discovery of el ectron wave filtering in Fe(100)/MgO 
crystalline systems^ 21141 where TMR up to 604% for a 
regular MT and 1056% for a double MT.pl were ob- 
tained. In the case of magnetic insulators, since the bar- 
rier height is spin-dependent, the current polarization in- 
creases with the barrier thickness^ P w tanhgd^. Here 
d, q, U and A are the barrier thickness, evanescent wave 
vector, barrier height and spin splitting of the magnetic 
insulator. Therefore, magnetic insulators have also been 
used as spin filters in quasi-spin valves structures. In 
these structures, the polarizer consists of ferro- or fer- 



rimagnetic tunnel barriers such as EuSp^jD NiFe20j2H, 
CoFe 2 Oj21 BiMnO^SI an d EuO 24 , the spacer is either 
tunneling (AlOx) or metallic and the analyzer is either 
a half-metal (LSMO), a transition metal ferromagnet or 
even another magnetic insulatoi^l. The magnetic insula- 
tors directly filter the tunneling electron spin and provide 
large TMR ratio at low temperature, up to 60% at IK in 
the case of fully insulating MTJsP. 

The mechanisms of spin-dependent tunneling through 
such structures have been investigated experimentally 
and theoretically unraveling complex bias dependence 
of the spin transport such as a spin- select ive Fowler- 
NordheirrPS (FN) tunneling processe d 18 ' 25 ' At large 
enough bias, the majority spin sees a reduced barrier 
thickness (FN tunneling) while the minority spin are still 
in the direct tunneling regime. The onset of the spin- 
selective FN process results in an effective enhancement 
of the TMH 2 ^, in sharp contrast with the monotonous 
decrease observed in reg ular MTJsP. 

In the context of the recent observation of spin transfer 
torque in YIGP^l, a few theoretical studies have addressed 
the nature of the spin torque exerted on magnetic insu- 
lators. Inoue^ derived the analytical expression of the 
electrically driven spin torque between a ferromagnetic 
metal and a magnetic insulator and Jia et alP3 calculated 
the mixing conductance at the interface between YIG and 
Ag using Density Function Theory. These two theoretical 
studies intend to explain the observed magnitude of the 
spin torque induced by spin Hall effect observed in Refs. 
[3] and 0] They show that although spin torque is small 
in such structures it can still be detected through spin 
wave measurements. In the present study, we focus our 
attention on the fully insulating MTJ proposed by Miao 
et alP. In particular, whereas Refs. [3J and E] measured 
the spin torque in the linear transport regime (low bias) , 
such fully insulating MTJ display strong non-linearities 
at large biases, illustrated by the dramatic enhancement 
of the TMR. This enhancement indicates a regime where 
the spin polarization is maximum and therefore one can 
reasonably expect a large spin transfer torque efficiency. 
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We show that the spin torque is dominated by the out- 
of-plane torque and can be experimentally exploited in 
current-driven magnetization dynamics studies. In Sec- 
tion [TTJ we discuss the nature of the spin torque in mag- 
netic insulators using a simple bilayer model. The theo- 
retical model for the fully insulating MTJ is exposed in 
Section [rn| and numerical results are reported on Section 
IV| Materials consideration and experimental relevance 
are discussed in Section |V| The conclusion is given in 
Section EH 



II. BILAYER MODEL 

Before entering in the details of the device physics, let 
us turn our attention to the spin transport taking place 
at the interface between two ferromagnets. This model 
follows the spirit of Stiles and Zangwill^ seminal work on 
spin torque in metallic spin-valves. Our system consists 
of two adjacent semi-infinite (either metallic or insulat- 
ing) ferromagnets whose magnetizations are misaligned 
by an angle 9. Considering a majority spin-polarized 
electron (f) flowing from left to right, its free-electron 
wave function reads 



il>(z < 0) = {e< z + R t 



-ik\ 



ip(z > 0) = r t e**«*| t) + Tje'***! |) 
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) — cos |||). iZf,4. is the spin-dependent reflection 

(transmission) coefficient and kj are the spin-dependent 
wavevectors in the ith layer. Note that these wave- 
functions are expressed in the quantization axis of the 
right layer. The solution of the system is obtained by a 
standard wave matching procedure throughout the junc- 
tion. The applied bias voltage V drives a (conserved) 
charge current J e and (spatially dependent) spin current 
J a through the junction. The charge and spin current 
densities carried by one electron described by the wave 
function ip are given by 
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where cr = {a x , a y , <x z ) is the vector of Pauli spin ma- 
trices. In the absence of spin-flip scattering, since the 
total angular momentum is conserved, the spin current 
lost is transferred to the local magnetic moment: the 
spin torque is directly related to the transverse spin cur- 
rent components J\ \ -\-iJ\_ giving rise to the in-plane and 
perpendicular torques, respectively. In the case of an in- 
terface between two metallic ferromagnets, one obtains 



the transverse spin current impinging on the right layer 
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To obtain the actual spin torque exerted on the right 
layer, the above expressions have to be averaged over 
the Fermi surface. In the case of metallic interfaces, it 
turns out that the perpendicular torque vanishes (a de- 
tailed discussion can be found in Ref. [SJ). At this stage, 
we point out that, due to the coherent interference of 
precessing spins, the torque extends in the right ferro- 
magnet on the distance A v w 2tt /(k R — k R f^. When the 
right layer becomes half-metallic, minority electrons can 
not propagate and an exponentially damped oscillating 
torque is founds 

The situation is quite different when considering the 
interface between a metallic ferromagnet and a insulating 
ferromagnet. In this case the wavevectors in Eq. [2] should 

The spin current reads 



be replaced by k R ^ 
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Since there are only evanescent states in the barrier, the 
torque exponentially decreases away from the interface. 
Interestingly, the perpendicular torque is now propor- 
tional to q R — q R which indicates that it should be smaller 
than the in-plane torque. Finally, in the case of two in- 
sulators, we obtain 
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(9) 
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Interestingly, there is simply no in-plane torque in this 
system and only a field-like perpendicular torque sur- 
vives. This can be understood by noticing that in the 
absence of propagating states, no spin precession occurs 
in the system. Therefore, the right layer feels an effective 
evanescent spin density aligned along z, which acts as an 
effective field. The more realistic case of a fully insulat- 
ing MTJ, comprising two magnetic insulators separated 
by a normal insulator, is more complex. However, the 
present simple interfacial models indicate that a domi- 
nant field-like torque can be reasonably expected in such 
systems. 



III. FULL JUNCTION MODEL 

We consider a fully insulating magnetic tunnel junc- 
tion as sketched in Fig. [l] in which two normal metals 



3 




FIG. 1. Schematics of the fully insulating MTJ. NM, MI, 
and S represent the non-magnetic electrodes, the magnetic 
insulators, and the insulating spacer, respectively. 
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TABLE I. Curie temperature (T c ), energy gap (U) and ex- 
change splitting (2A) of several ferro- and ferrimagnetic insu- 
lators. 



(NM) work as electrodes and two magnetic insulating 
layers (MI) are separated by one non-magnetic insulat- 
ing spacer (S). The magnetization of the left magnetic 
insulator Mi is aligned along z and the one of the right 
magnetic insulator M R lies in the x-z plane and is rotated 
by an angle 9 with respect to Mi. 



Magnetic insulators are present under various crys- 
tal structures such as rocksalt (EuS, EuO), perovskite 
(BiMnOa), inverse spinel (NiFe204,CoFe204) or garnet 
(Y 3 Fe50i2) and their Curie temperature varies over a 
wide range (16.6-850K). They display a gap of the or- 
der of U «1±0.5 eV and their conduction band shows 
spin splitting of the order of A w 0.3 — 2eV (see Ta- 
ble [I]) . In the present model, the band structure of these 
materials is represented by a simple free electron model. 
Although important features such as interfacial resonant 
states and electron wave filtering might have some im- 
pact on the tunneling spin transport (as in MgO/Fe(100) 
systems^!), we disregard these effects at this stage and 
retain only a minimal set of parameters. The barriers 
height, thickness and splitting of the insulators are de- 
noted \Jl,s,r> dh,s,R, and ^l,r respectively. Note that 
the spacer is a non-magnetic insulator (As = 0). Under 
the free electron approximation, the Hamiltonian of the 
system reads 



H 



d 2 



2m* dx 2 
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where d is the total thickness of the composite barrier, to* 
is effective mass, and Ep is Fermi energy in the metallic 



electrode. The solution of Eq. (11) is on the form 



|V) = Y} aA i(Z°ix)) + bB^ixMa) (12) 



where Ai and Bi are the homogeneous Airy functions, 
and the argument Z a (x) reads 



dh 2 



to* I eV | 



2/3 
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where E, k and a — ± are the energy, the in-plane 
wave vector and spin projection of the incident electron. 
Previous methods such as Refs. US] and have used 
WentzclKramersBrillouin (WKB) approximation instead 
of the exact Airy functions to solve the eigenstates of Eq. 
(111. The WKB approximation is justified in the case of 



thick barriers only and does not capture the details of 
the Fowler-Nordheim resonances when the applied bias 
voltage exceeds the barrier height. On the other hand, 
whereas Airy functions are the exact solutions of the 
free electron Hamiltonian, they are numerically stable 
for only thin enough barriers. Thus, the gain in accuracy 
in describing the physical tunneling process at large bi- 
ases is balanced by the constraint on the barrier thickness 
and height. Therefore, in the present study we choose to 
use the exact solutions provided by Airy functions but 
to apply this method to thin barriers only. Although the 
thickness of the barriers studied are smaller than the one 
usually attainable in the experiments, we believe that 
the physics addressed in this article is better described 
by the present method and will not significantly change 
for thicker barriers. 

The spin-polarized current flowing throug h the struc- 
ture exerts a local spin transfer torque (STTp^ on the 
magnetic layers. The torque estimated in the left MI 
layer per unit area is equal to 

T=-Jv- J s = J s {d L )- J 3 (0) (14) 
= T n M L x (M L x M R ) + Tj_Ml x M R (15) 

where Tj | is the in-plane torque and T± is the out-of-plane 
torque. 



IV. RESULTS 

Unless stated otherwise, the parameters for the system 
are: U l ,r= 0.8 eV, A l ,r= 0.1 eV and U s = 1 eV. In our 
definition, at positive bias, electrons travel from the left 
electrode to the right electrode. 
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FIG. 2. Left panel: Bias dependence of the TMR for (a) 
different spacer thicknesses (b) different spacer barrier height 
and (c) different exchange splitting. Right panel: Energy 
profile of the fully insulating magnetic tunnel junction: (d) A 
direct tunneling occurs when the applied voltage is less than 
the lowest spin-up barrier height; (e) The Fowler-Nordheim 
tunneling is observed when the bias voltage exceeds the lowest 
spin-up barrier heiglrP^. 



A. Tunneling Magnetoresistance 

The bias dependence of the TMR, defined as as 
TMR=(J P -J A p)/ Jap (where P and AP refer to the par- 
allel and antiparallel states, respectively), is represented 
in Fig.2 (a)-(c). As observed by Miao et alPJ, the TMR 
displays peaks at both positive and negative bias polari- 
ties. The feature has been identified by Miao et aiP^ as 
a direct consequence of so-called Fowler-Nordheim res- 
onant tunneling^!. When the applied bias voltage ex- 
ceeds the barrier height for the majority spins, the ef- 
fective thickness of the potential barrier seen by these 
electrons is reduced whereas it is unchanged for the mi- 
nority spins. Therefore, the current for majority spins 
dramatically increases resulting in a significant enhance- 
ment of the TMR [see Fig.2(d,e)]. Pig.2(c) illustrates 
this effect for different exchange splitting A. As ex- 
pected, the effect is stronger for larger exchange splitting, 
due to a stronger spin filtering. Similarly, increasing the 
spacer thickness [Fig. 2(a)] or decreasing the spacer bar- 
rier height [Fig. 2(b)] enhances the TMR peak which oc- 
curs at smaller bias. 

Interestingly, at larger bias a reversal of the TMR is ob- 
served as shown in Fig. 6(a). When further increasing the 
bias voltage, the TMR ratio exhibits a clear oscillatory 



behavior and the period of the oscillations increase with 
the bias voltage. This oscillation has also been obtained 
by Montaigne et alP^ is a regular MTJ when using an 
exact Airy description of the tunneling transport. This 
is directly related to the resonances associated with the 
FN process. Note that t hese o scillations are not captured 
by WKB approximation^^!. 



B. Spin Transfer Torque 

To investigate the nature of the spin transfer torque, 
we set the magnetizations in the perpendicular configu- 
ration, (6 = 7r/2). The bias dependence at low voltage 
for both components of the spin torque is presented in 
Fig. 3. The form of the spin torque is 



7|| = aiV + a 2 V 2 + a 3 V 3 
T± = b + b 1 V + b 2 V 2 



(16) 
(17) 



The in-plane torque is essentially antisymmetric in bias 
voltage (note a small asymmetry at large bias) whereas 
the out-of-plane torque is essentially symmetric present- 
ing a quadratic bias dependence (with, again, a small de- 
viation at large bias). This is in sharp contrast with the 
bias dependence of the torque obtained in regular sym- 
metric magnetic tunnel junctiomP^^. For comparison, 
Theodonis et alW^ found a dependence on the form 



7|| 



= aiV + a 2 V 2 
= b + b 2 V 2 



(18) 
(19) 



Furthermore, whereas the out-of-plane torque is found to 
be smaller than the in-plane torque in conventional tun- 
nel junctions^, we find that it dominates the spin torque 
in fully insulating junctions. The influence of the bar- 
rier thickness and height are reported in Fig. 3 (a,c) and 
(b,d), respectively. The magnitude of both in-plane and 
out-plane torques decrease when increasing the spacer 
thickness or height, as expected. The torque magnitude 
as a function of the exchange are plotted in Fig. 3 (e,f). 

At large biases, the spin torque shows an dramatic in- 
crease by 2 to 3 orders of magnitude, as displayed in Fig. 
4. Most importantly, whereas the out-of-plane torque is 
smaller than the in-plane torque at small biases (see Fig. 
3), it largely exceeds the in-plane torque at large biases. 
The contribution of majority and minority electrons orig- 
inating from left and right reservoirs to the spin torque 
has been calculated in Fig. 5. Whereas the torque orig- 
inated from majority and minority spins compete with 
each other to produce the effective in-plane torque, they 
add up to produce the out-of-plane torque. In particu- 
lar, note that the resulting spin torque is more efficient 
at negative bias, which corresponds to a direct torque, 
the electrons flowing from the right electrode to the left 
one. 

Mirroring the TMR behavior, the spin transfer torque 
at very large bias displays an oscillating bias dependence 
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FIG. 3. Low bias dependence of (a,d) Tj| and (b,d) T± com- 
ponents of the torque exerted on left MI layer in a symmetric 
junction for 9 — n/2 together with the corresponding effects 
of different spacer thicknesses (a,c), different spacer barrier 
heights (b,d) and different exchange splitting (e,f). 



as shown in Fig.6(b,ch somewhat similar to the one ob- 
served by Tang et air** in regular MT Js. This oscillation 
is attributed to the resonant states appearing when the 
barrier height goes below the Fermi energy of the incom- 
ing electrons, as mentioned above. 



V. DISCUSSION 

The magnitude of the spin transfer torque calculated 
using the present model is w0.5 eV.[im~ 2 at IV, which 
is about two orders of magnitude smaller than what has 
been calculated in regular MTJs using the same free elec- 
tron approach (30-100 eV [imT 2 at IV in Ref. HO]). How- 
ever, in the present discussion we show that in principle, 
this difference does not prevent from observing a measur- 
able effect of the current on the magnetization dynamics. 

First, we have deliberately chosen conservative mate- 
rials parameters, such as a small exchange energy A = 
O.leV. As shown in Fig. 3(e), the torque is more-than- 
linear as a function of the exchange and therefore, in real- 
istic materials a larger spin torque is expected (see Table 
I). From the materials viewpoint, it is necessary to search 
for compounds with low barrier height, strong spin split- 
ting and large Curie temperature (to favor thermal acti- 
vation of the magnetization dynamics). Following Table 
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FIG. 4. High bias dependence of (a,b,c) T\\ and (d,e,f) T±, to- 
gether with the corresponding effects of (a,d) different spacer 
thicknesses (b,e) different spacer barrier heights (c,f) and dif- 
ferent exchange splitting. 



|TJ CoFe 2 04 might be an interesting candidate. However, 
a hurdle remains the thin layer constraint. Miao et al.^ 
were able to grow 0.6 nm thick of EuS layer by thermal 
evaporation, which renders spin transfer studies reach- 
able. Furthermore, the AlOx spacer could be replaced 
by a metallic spacer. In such a structure, resonant states 
would develop in the metallic spacer that could enhance 
the torque efficiency 

The fully insulating MTJ structure studied in this work 
is relevant for two types of experimental configuration. 
The first one is the usual bias voltage-driven magneti- 
zation switching such as the one currently observed in 
regular MTJa^M In this configuration, when the magne- 
tization of the layers lies in the plane, the intrinsic crit- 
ical switching voltagefat zero temperature) is denned 
V c = ±^(Hk + Hdj^i, where d is the layer thickness, 
Hk and Hd are the anisotropy and demagnetizing fields 
and cii is the spin torque efficiency, defined in Eq. (16 1. 



Another configuration is the bias voltage-driven domain 
wall motion, as proposed in Ref. 031 In this case, the 
spin transfer torque is used to move a domain wall con- 
tained in the free layer. Whereas the domain wall de- 
pinning is governed by the in-plane torque, the domain 
wall steady state velocity is governed by the out-of-plane 
torque v s» — ^A. The interesting feature is that both 
quantities, magnetization switching threshold and do- 
main wall velocity, are governed by the ratio between the 
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FIG. 5. Spin-dependent and reservoir-dependent contribu- 
tions to the in-plane (a) and out-of-plane torques (b), exerted 
on the left MI layer. 
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It is however important to notice that the extremely 
small magnetic damping of YIG is observed only for thick 
layers (>60 nm). Decreasing the layer thickness results 
in an increase of interfacial spin pumping effectP (and 
therefore an increase of the effective damping), as well 
as a modification of the band structure of the magnetic 
insulator^ (associated with a decrease of the Curie tem- 
perature) . Finally, YIG is known for its vanishing damp- 
ing and stronger ferromagnets such as EuO or CoFe2C"4 
have generally larger damping factors. This may result 
in weaker current-driven magnetization dynamics which 
can still be detected using methods such as spin diode 
and spin-FMR measurements^. 



VI. CONCLUSION 

The current-driven spin torque in fully insulating tun- 
nel junctions composed of ferromagnetic insulating elec- 
trodes has been studied theoretically using a free electron 
description. We find that: (i) the out-of-plane torque 
and in-plane torque are symmetric and antisymmetric at 
low bias voltage, respectively (ii) both of torques show a 
dramatic enhancement at large bias voltage (iii) the out- 
of-plane torque dominates in-plane torque at large biases, 
(iv) at very high bias voltages, the TMR and spin torque 
shows an oscillatory dependence. Finally, we discussed 
the experimental conditions of the observation of current- 
driven magnetic domain wall motion and magnetization 
switching. 
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FIG. 6. Bias dependence of the TMR (a), in- plane (b) and 
out-of-plane torques (a) at very large biases for different 
spacer thicknesses. 



spin torque and the damping constant a. As mentioned 
in the introduction, since electron-magnon interactions 
are quenched in magnetic insulators, the damping is usu- 
ally about much smaller in magnetic insulators than in 
transition metals (a ss 6.7 x 10~ 5 for YIG compared to 
0.01 for CoFe). Consequently, one can expect that the 
efficiency of the spin torque needs not to be as large as in 
MgO-based MTJs to observe the current-driven magneti- 
zation dynamics. Following our crude estimations, a bias 
voltage of 1 V should be sufficient to observe the magne- 
tization dynamics in the case of a damping of a « 10~ 4 . 
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